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CONTINUOUS INVENTORY 

"Today I have grovn taller from walking 
with the trees." 

...Karle Wilson 


Milwaukee, Wls. January, 1964 Mo. 118 


The New Year, the New Year. Everywhere 
the New Year. The Old Year was already 
looked upon as dead; and its effects were 
selling cheap, like some dead mariner* s 
aboard ship o Its patterns were Last Year* s 
and going at a sacrifice, before its 
breath was gone. Its treasures were mere 
dirt, beside the riches of its unborn 
successor. 


CHARLES JOHN HUFFHAM DICKENS 
The Chimes. Fourth Quarter. 
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A Statistical Check Computer Program for the I»B.M. 1620 


C.F.I. can provide the forester and business manager with a vast number 
of results. Every inventory is only a sample of the total forest, and 
results obtained from inventory data analyses must be inferred from the 
sample taken. The results are only approximations, however, but it is 
assumed when taking the sample that data gathered are representative 
of the total forest. Consequently, a risk of being wrong is accepted 
when results are stated. What is required as a supplement to the results 
is a measure of their reliability. The Statistical Check provides a 
means to determine this reliability. Our Newsletter this month deals 
with a Statistical Check procedure that should be used on all inventories 
and it is the first in a series of Newsletters deali n g with computer 
applications to the Region 9 C.P.I. Job. 

The computer program presented was written to be used on I.B.M. 1620 
computers. The 1620 was chosen particularly because many universities 
possess installations of this type which can he used by those operating 
a C.P.I. system. This program is not, however, restricted to a 1620. 

With a few minor modifications, which will be apparent to a computer 
programmer, the program, written for FORTRAN II, can be compiled and 
run on larger computers such as the I.B.M. 7090* Prepared FORTRAN 
source decks of the program are available upon request at the U. S. 

Forest Service Regional Office In Milwaukee. 

The compiled program consists of one object deck which will compute 
area break statistical checks on volume (or growth) estimates, area 
estimates and volume and area estimates combined, and punch results. 

For various area breaks common to C.F.I. in Region 9 be** 1 pulpwood 
and sawlog quantitative input data can be used. Fourteen area breaks 
segregated into administrative, operational, physical-graphical and 
ecological groups can also be included in the input data. A separate 
run is required for pulpwood or sawlog input data, although it is not 
necessary to reload the program ctbject deck after each run. When an 
analysis is completed the progr am brandies to the beginning of the 
routine and cells for more input data to be analysed. 

Specification of pulpwood or sawlog data is made by using a Parameter 
Card which precedes the input data when the latter are read into the 
machine. The Parameter Card layout appears on page 5 » and instructions 
for its use are presented with the program operating instructions. 

Input and output card layouts for pulpwood and sawlog data also appear 
on this page. The only input card field titles that need conform with 
those that appear In the card layouts ares (l) the number of plots 
in a particular area break; (2) the sum of the plot volumes; 

(3) the sum of the plot volumes squared, and (k) the expanded area. 

It is suggested that the area break columns for input cards be used 
as designated on the card layouts; there is a fixed correspondence 
between input card columns and output card columns. 
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The original object deck was compiled using FORTRAN IX on a 60,000 
bit storage 1620 with the Indirect addressing feature and automatic 
floating point multiply and divide. A listing of the source deck 
from which program object decks are generated appears on page 6, 
and a corresponding flow diagram is presented on page T« 

Statistical Procedure Leaflet #9 which accompanied C.F.I. Newsletter 
#79> October, i 960 contains the formulas used in the Statistical 
Check Program. All results are calculated at a confidence level 
of approximately 95$ with the assumption that sample plots are 
l/5-acre in size. The numerical constants in source card numbers 
00037 and 00075 for other plot sizes should be changed, for example, 
frcm 5.0 for a l/5-acre plot to 7.0 for a l/7-acre plot. 

Limits of error are computed as non-percentage figures and must be 
multiplied by 100 to yield percentage results. When the computed 
results exceed the specified size of the output fields, the 
machine types and punches an error message. 


James N. Hool 

Department of Forestry & Conservation 
Purdue University 
Lafayette, Indiana 


This program was first used under 
the direction of the consulting 
firm of George Banzhaf & Company, 
Milwaukee, Wisconsin. Hie CFI 
case involved 1,100 pexmanent 
plots of fixed radius, establish- 
ed on the timber lands of the 
Cleveland-Cliffs Iron Co. in 19^1 



rating Instructions 



3. Load program object deck: 


4. Load program subroutine deck: 
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Machine Operating Instructions (continued) 

5« Enter data to be processed: Place Parameter Card aheefl of 

data Input cards. 

Place Parameter Card and data 
Input cards in reader feed. 

Press START on console 

Press READER START and PUNCH START 

Press READER START to read 
last data input card. 

6. Remove output data from punch hopper and list. 
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STATISTICAL CHECK FORERAN SOURCE-DECK LISTIHG 


6 


Card Columns 
Identification 


Initialization 


Computation 

end 

Output 


I 


Computation 

and 

Output 


If number of 
plots Is less 
than f lye 

Computation 

If £ plot 
volume equals 


COMBINED PULPWOOD-SAWLOG STATISTICAL CHECK ROUTINE FOR FOREST I NV 
C.F.I. DIVISION, S-PF, REGION 9, USF S-M I L w AUKEE * W I S 


I PLOT S = 3 

10 READ 20 » 1 L 1 * FN, 1L2 

20 FORMAT (II, IX, F4.0, IX, 14) 
IF ( IL1 - 1 ) 30, 30, 210 


30 PUNCH AO 

AO FORMAT ( 2 7H.STAT I ST I C AL CHECK* PULPWOOD) 

I = 0 

110 READ 60, I A 1 , 101* 102, IP1, IP2, IP3, IPA, IE1, IE2, IE3, 

DIES, I E 6 , I E 7 , I C 1 , IC2, NA, SX , SOX, EXPA 
60 FORMAT (611, 13, 12, 311, 212, 211, 12, I A', IX, F9.3, IX, 

1F9.3, IX, F 7 • 0 ) 

IF (NA - 3) 70, 120, 120 

70 ~ PUNCH 80, I A 1 » 101, 102, IP1, IP2, IP3, IPA, IE1, IE2, IE3, 
1 1 E 5 * I E 6 , I E 7 , IC1, IC2, NA, EXPA 
80 FORMAT (II* 2X, II. IX, II, 2X. II. IX, II, IX, II, IX, 13, 


1 I E 5 , I E 6 , I E 7 , IC1, IC2, NA, EXPA 
80 FORMAT (II. 2X, II, IX, II, 2X, II, IX , I 
112, IX, II, IX, II, IX, II, 12* 12, II, 2 
227H NO. OF PLOTS LESS THAN 5, 3X, F9.0) 
1 = 1 + 1 

IF ( IL2 - I ) 390, 390, 110 
120 FN A = NA 

T 1 = ( FN - FN A ) / ( FN * FN A ) 

EA = 2.0 * SQRTF(Tl) 

IF (SX) 170* 130* 170 

TTa'PUNCH 1 A 0 , i A 1 , 101, 102, TpTT IP2, IP3, 


11, IX, II, IX, II, IX, I 3 1 

12 , II, 2 X , II, 12 , IA, 




I E 3 , 

i I E 6 

, 

It/, 

1C1 , 

IC2 

FORMAT ( 

I 1 

, 2 X , 

- 11, 

IX , 

12 , 

IX, 

I 1 

* 1 X j 

. 11, 

IX , 

11H 

VOL 

• 

ZERO, 

> F6.3, 1 


, IP1, IP2, I P 3 , IPA, I E 1 , IE2, I E 3 
NA, EA, EXPA 

II, 2 X , II, IX, II, IX, 11, IX, 13, 
II, 12, 12, II, 2 X , II, 12, IA, 


1 = 1 + 1 

IF ( IL2 ~ I) 390, 390, 1 10 

1 70 ~ B = FNA * SOX/SX**2 - 1.0 

EV = 2.0 * SQRTF ( B/ ( FNA - 1.0)) 

CSQ = B * FNA/ ( FNA - 1.0) 

EVA = 2.0 * SQRTF ( B/ ( FNA - 1.0) + Tl) 

VPERA = 3.0 * SX/FNA 

PUNCH 180, I A 1 * 101, 102, IP1, IP2, IP3, IPA, I E 1 . IE2, IE3 

1 I E 3 , IL6, I E 7 , IC1, IC2, NA, CSQ, LV , EA, EVA, VPERA, EXPA 
180 FORMAT (II, 2X, II, IX, II, 2X, 11, IX, II, IX, II, IX, 13, 
112, IX, II, IX, II, IX, II, 12, 12, II, 2 X , II, 12, IA, 
2F3.2* F6.3, K6.3, F6.3, F7.3, F9.0) 

1 = 1 + 1 

IF ( IL2 - I ) 390, 390, 110 



imsiiii 


If number of 
plots Is less 
than five 

Computation 

If € plot 
volume equals 


II, IX, II, IX, 13, 2 X , 


FORMAT ( 2 5HSTAT I ST I CAL CHECK, SAWLOG ) 

I = 0 

READ 2 AO , I A 1 , 101, 102, I P 1 , IP2, IP3, IP4, I El , IE2, IE3, I E* 
1IE5, 1E6, I E 7 , IC1, IC2, NA, SX * SOX, EXPA 
FORMAT (611, 13, 12, 311, 212, 211, 12, IA, IX, F8.1, IX, 

1 F 1 1 . 1 , IX, F 7 . 0 ) 

IF (NA ~ 5) 230, 300 , 300 

"PUNCH 260, I A 1 , 101, 102, I P 1 , IP2, IP3, IPA, IE1, IE2, IE3» IE 
1 I E 5 , IE6, I E 7 , I C 1 * IC2, NA, EXPA 
FORMAT (II, 2 X , II, IX, II, 2X, II, IX, II, IX, II, IX, 13, 2 X , 

112, IX, II, IX, II*. IX, II. 12, 12, I 1 , . 2 X , II, 12, IA, 

e. 2 7 H NO. OF PLOTS LESS THAN 3, 3X, F9.0) 

1 = 1 + 1 

IF ( IL2 - I ) 390, 390 , 290 
"FNA = NA " ‘ 

Tl = ( FN - F NA ) / ( F N * FNA) 

EA = 2.0 * SORTF ( T 1 ) 

IF (SX) 3 30 , 3 10, 3 30 

"PUNCH 320 , I A 1 , 101, 10^* IP1, IP2, IP3* IPA, IE1, IE2* IE3, lb 
1 I E 3 , I E 6 , 1 E 7 , I C 1 , IC2, NA, EA, EXPA 

FORMAT (II, 2 X , II, IX, II, 2X, II, IX, II, IX, II, IX, 13, 2X, 

112, IX, II, IX, II, IX, II, 12, 12, 11, 2 X , II, 12, IA, 

niH VOL. ZERO, F6.3, 13X, F9.0) 

1 = 1 + 1 . 

IF ( IL2 ~ I ) 390, 390, 290 

"B = FNA * SOX/SX#*2 - 1.0 
EV = 2.0 * SOR T F ( B / ( EN A - 1.0)) 

CSQ = B * FNA/ (FNA - 1.0) 

EVA = 2.0 * SQRTF (B/( ENA - 1.0) + Tl) 

VPERA = 3.0 * SX/FNA 

PUNCH 360, I A 1 , 101, 10^, I P 1 , IP2, IP3, IPA, IE1, IE2, IE3, U 


IX, II, IX, 111 
II, 2 X , II, 12 , 


IE2, 

IE3- 

> IE 

IX , 

13, 

2 X , 

I A, 




CSQ, EV, EA, EVA 


I E 1 , IE2, I E 3 , 
■ VPERA, EXPA 



1 I L 3 , I E 6 , I L 7 , IC1, I C 2 , NA, CSQ, EV, EA, EVA, VPERA, EXPA 

FORMAT (11, 2 X , II, IX, II, 2X, II, IX, 11, IX, II, IX, 13, 2X - 

112, IX, II, IX, II, IX, II, 12, 12, II, 2 X , II, 12, IA, 

^ F 3 • 2 , F 6 • 3 , F6.3, F 6 • 3 , h7,0, F 9 • 0 ) 

1 = 1 + 1 

IF I IL2 - I ) 390, 390 , 290 

"TYPE ~ ~ 

FORMAT ( 3 1 H ENTER PARAMETER CARD AND DATA PRESS READER START) 

GO TO 10 









JAMES M. WOOL 
12 - 13 - 63 











FLOW DIGRAM FOR 


\$ZO STATISTICAL CHECK PROGRAM 



BLOCK* SOURCE CARD* BLOCK ; 
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STATISTICAL PROCEDURE LEAFLET #9 
CALCULATING THE LIMIT OF ERROR IN AREA AND VOLUME COMBINED 


Since there is a limit of error for both the estimated average volume and the 
estimated area, when we multiply the two to get the total volume on any 
particular breakdown we shall wish to calculate the limit of error for the 
product. The limit of error for this product is very easily calculated by 
adding data used to calculate the individual limits of error in volume and 
area as shown below. 

Formulae for limits of error 1 /: 


In VOLUME 


In AREA 


In combined AREA x VOLUME 


Ev 


i2 Vi 


B 


Na - 1 


go _ + 

“ 2 V E x Na 


' i2 VsrVf 


- Na 


x Na 


WHERE: 


Bv 

Ba 

Eva 

N 

Na 

Ex 

Ex 2 

(Zx ) 2 

B 


For EXAMPLE when: 


Limit of error-VOLUME 


Ev = + 2 * / ■§ 10 
V 31-1 


limit of error in volume 
limit of error in area 

limit of error in combined area and volume 

total number of plots on whole forest 

number of plots in an area breakdown 

sum of individual plot volumes 

sum of squares of individual plot volumes 

square of the sum of individual plot volumes 


Na (gx ) 


(Ex)‘ 


N 

Na 

B 


66l plots 
31 plots 

31 (93.081) 
1594.002 


rx 2 

(rx ) 2 


* 93-081 

= 1594.002 


Limit of error -AREA 


Ea = + 2 -\ f 

- V 661x31 


1, or .810 
Limit of error-AREA x VOLUME 


^ ^ _ r&io . 661-31 
** * ±2 \ 31-1 + 66S31 


12 V" 


810 

30 


± 2 y. 027000 

+ 2 (.164) 


1 2 


}0_ 

£1 


+ 2 -\/. 030745^ 
± 2 (.175) 


f .810 . -636" 
- 2 V 30 + 20I91 

+ 2 y.0577^5 

+ 2 (.240) 


± -328 , or ± 33* 


± .350 , or ± 35* 


+ .480, or + 48* 


1/ The 2 in each formula is approximate''^ factor for 95* probability. See 
Statistical Procedures Leaflet #4 for discussion of probability. 

WILLIAM WARREN BARTON, Forester,, U.S. Forest Service, Region 7 





